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Oxide heterostructures with polar films display special electronic properties, such as the electronic 
reconstruction at their internal interfaces with the formation of two-dimensional metallic states. 
Moreover, the electrical field from the polar layers is inversion-symmetry breaking and generates a 
Rashba spin-orbit coupling (RSOC) in the interfacial electronic system. We investigate the quantum 
capacitance of a heterostructure in which a sizeable RSOC at a metallic interface is controlled by 
the electric field of a surface electrode. Such a structure is, for example, given by a LaAlOa film on 
a SrTiOa substrate which is gated by a top electrode. Such heterostructures can exhibit a strong 
enhancement of their capacitance [I]. The capacitance is related to the electronic compressibility 
of the heterostructure, but the two quantities are not equivalent. In fact, the transfer of charge to 
the interface controls the relation between capacitance and compressibility. We find that due to a 
strong RSOC, the quantum capacitance can be larger than the classical geometric value. However, 
in contrast to the results of recent investigations [9llll| the compressibility does not become negative 
for realistic parameter values for LaAlOs/SrTiOs and, therefore, we find that no phase-separated 
state is induced by the strong RSOC at these interfaces. 

PACS numbers: 71.70.Ej, 73.20.-r, 73.40.Rw 


I. INTRODUCTION 

The static capacitance of a complex metallic structure 
is commonly considered to be a global quantity of the en¬ 
tire system and, consequently, is supposed to characterize 
the electronic state of the system rather insufficiently to 
deduce microscopic properties. However this is not al¬ 
ways the case. For instance, recent experiments on polar 
oxide heterostructures with conducting LaAlOa/SrTiOa 
interfaces suggest that the capacitance for these systems 
is a quantity which allows to draw distinct conclusions on 
the electronic interface state mm- A special character¬ 
istic of this interface is that the electrons are subject to 
a sizeable spin-orbit coupling [Si]. The question arises 
if spin-orbit coupling is a relevant control parameter for 
the capacitance of the interface. 

Here, we address the fundamental issue whether spin- 
orbit coupling of the interface electrons allows to mod¬ 
ify the capacitance of the heterostructure significantly, 
and we estimate the spin-orbit contribution for systems 
based on LaAlOa/SrTiOa interfaces. An enhancement 
of the capacitance is possible in the presence of a suf¬ 
ficiently large Rashba-like spin-orbit coupling (RSCO), 
because the RSOC allows to alter the electronic band 
energies by controlling the charge carrier density. Such 
a RSOC appears to explain magnetotransport measure¬ 
ments at LaAlOa/SrTiOa interfaces [31I11|BHH]) and it has 
been shown within ab initio evaluations that a RSOC is 
generated by the electric held of the polar LaAlOa layers 
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on top of the SrTiOa substrate (5). The RSOC lowers 
the band edge quadratically with the spin orbit coupling 
parameter or (see Fig. [^. This even allows for a neg¬ 
ative compressibility assuming that or depends on the 
electric held E of the polar LaAlOa layers which in turn 
is related to the charge density at the LaAlOa/SrTiOa 
interface mm- With their spectacular result of a nega¬ 
tive compressibility in the presence of a sizeable RSOC, 
the authors of Refs. |9HII] suggest that this provides 
an intrinisic mechanism for the inhomogeneous phases 
observed at LaAlOa/SrTiOa interfaces—yet, other sce¬ 
narios for phase separation at LaAlOa/SrTiOs interfaces 
have also been investigated recently mm- 

It is well known, that a negative electronic compress¬ 
ibility on one or both electrodes of a capacitor can lead 
to an enhancement of the capacitance beyond its geo¬ 
metrical value m [IS]. In fact, the geometric capaci¬ 
tance of a two-plate capacitor, Cgeom = sA/Aird, is mod- 
ihed by electronic contributions which are taken into 
account in the electronic compressibility of the plates 
Ki = (nf dfj,i/dni)~^, the derivative of the charge carrier 
density Ui at plate i with respect to the chemical poten¬ 
tial ^i- Here e is the dielectric constant of the material 
between the two electrodes of area A, and d is the thick¬ 
ness of the dielectric. The kinetic contributions to the 
capacitance were coined by Luryi [TB] as quantum capac¬ 
itance. This notation may be extended to all electronic 
compressibility terms (see Ref. US]). 

Already in the early nineties, a capacitance enhance¬ 
ment of semiconductor heterostructures has been iden- 
tihed when the electronic compressibility of quantum 
Hall systems was investigated [14]. The enhancement 
of the semiconductor capacitance could be traced quan- 
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titatively to the presence of quantum exchange effects 
in the electrodes. The contribution from electronic cor¬ 
relations to the capacitance is not negligible for very 
dilute electron systems m and for electronic systems 
close to half-filling However, the charge density at 
LaAlOa/SrTiOa interfaces is in the intermediate range, 
neither as low as in the semiconductor systems m nor 
close to half-filling. Therefore, the more sizeable capac¬ 
itance enhancement of the oxide LaAlOa/SrTiOa het¬ 
erostructure [1] and its negative compressibility [2] can¬ 
not be explained exclusively by exchange and correlation 
effects and other mechanisms might become relevant. 

To determine the thermodynamic electronic compress¬ 
ibility, two alternative procedures are conceivable: (i) one 
can either consider the interface subsystem with an ex¬ 
ternal electric field which induces the RSOC. In order 
to identify the electronic compressibility of the interface 
electrons, one has to keep the electric field hxed in ther¬ 
modynamic derivatives. Or (ii), one evaluates the com¬ 
pressibility for the entire electronic system, for which the 
RSOC-generating field is an internal field. 

A scheme in the sense of (i) has been suggested by 
the authors of Ref. [3]. There the electric field was set 
as a function of the interface charge carrier density n 
and, then, the inverse compressibility was identified from 
the derivative of the chemical potential fi with respect 
to n, keeping any dependence else fixed. We disagree 
with this approach because there the electric field has 
to be kept hxed and therefore does not contribute to the 
compressibility [18]. In this paper we establish the second 



FIG. 1. (Color online) Top: Dispersion of a planar electron 
system with RSOC. The dashed line indicates the dispersion 
for zero Rashba coupling. Bottom: Occupation in fc-space for 
fj .1 < 0 (left, only the lowest band is occupied) and for p 2 > 0 
(right, both bands are occupied). 


procedure (ii) and relate the electronic compressibility of 
the entire electronic system directly to the capacitance 
of the device—consisting of the interface, dielectrics, and 
the surface of the LaAlOa him or a metallic top electrode. 

As indicated above, previous investigations focused on 
the interface subsystem mm, whereas we consider the 
entire system. For realistic values of the RSOC, we iden¬ 
tify a positive compressibility for the most part of the 
parameter space indicating that the system is thermo¬ 
dynamically stable. Nevertheless, we hnd an increased 
capacitance with respect to its geometrical value, the ca¬ 
pacitance enhancement being caused by the RSOC in the 
metallic interface. 

In this article, we will investigate two different models, 
both building on a two-dimensional (2D) metallic system 
in which a RSOC modihes the electronic state: First we 
analyze a capacitor with two metallic planes, where the 
top electrode is a 2D or 3D metal (gate) and the second 
plane is the interface between a polar and non-polar insu¬ 
lator. The interface attains a metallic state with RSOC 
upon charge transfer from the top plate to the interface 
and the concomitant build-up of an internal electric field. 
Finally we study a two-layer system with the spin-orbit 
coupling controlled by the filling of three t 2 g-orbitals at 
the interface. We will show that the multilayered system 
has a range of parameters where it has a positive total 
compressibility, but its capacitance is enhanced above its 
geometrical value on account of the spin-orbit coupling 
in the metallic interface plane. In this analysis we ad¬ 
just the parameters to those specihc for LaAlOa/SrTiOa 
(LAO/STO) heterostructures in order to explore the fea¬ 
sibility of spin-orbit controlled capacitance enhancement 
with a polar dielectric material. 

II. RASHBA SPIN ORBIT COUPLING AND 
COMPRESSIBILITY 

In a 2D free electron system subject to a Rashba spin- 
orbit coupling (RSOC) 7(E x p) • S, with coupling con¬ 
stant 7 , the dispersion (k) is given by (see, for example. 

Ref. [ 11 ]): 

where Nq = mo/TrlT? = (mQ/me)/{Tra^e'^) is the density 
of states for both spin directions in the absence of RSOC, 
an defines the strength of the RSOC, and mo is the ef¬ 
fective mass of the charge carriers. The RSOC is made 
possible by the absence of inversion symmetry, caused 
by an electric field perpendicular to the two-dimensional 
layer containing the electrons with areal density n. The 
spin-orbit splitting of the dispersion, Eq. Q, is shown 
in Fig. [2 If /r < 0, only the lower band is occupied and 
the Fermi surface is an annulus; for /r > 0 the two occu¬ 
pied bands yield two concentric circles in fc-space. As the 
lowest energy = —ttNoo^ of the dispersion is pro¬ 
portional to the strength of the RSOC, it decreases with 
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increasing strength of the symmetry breaking electrical 
field. Thus a stronger electrical field lowers the chemi¬ 
cal potential for constant electron density n. In Ref. [5] 
it was suggested that a negative compressibility can be 
obtained from exploiting this effect. 

The inverse compressibility is obtained by the second 
derivative of the free energy F per area with respect to 
the electron density with all external variables held fixed. 
Since the symmetry breaking electrical field (SBEF) is an 
external variable in this scheme (and the field generated 
by the electrons in the plane themselves is symmetric), 
Or, is independent of n in the thermodynamic derivatives, 
and there is no dependence of the compressibility on the 
RSOC. 

Here we propose an extended scheme in which the 
electric field from the polar film, and consequently the 
RSOC, is an internal parameter, i.e., a quantity that de¬ 
pends on the electron density n. Then the density de¬ 
pendence of Or has to be accounted for when taking the 
derivative of the energy for the compressibility. In ad¬ 
dition to the interface layer with a finite RSOC (TO), 
a second layer (LI), e.g., the surface of the LAO film 
in the LAO/STO heterostructure, comprising electrons 
and a positive background charge, is placed parallel at a 
distance d. Eventually, both layers—each of them electri¬ 
cally neutral—are connected. Since typically both layers 
had different electrochemical potentials, electrons will be 
exchanged between LO and LI until both layers are on the 
same electrochemical potential. It is of no relevance for 
the static capacitance how electronic charge is transferred 
between the plates. However, an electric field, which is 
proportional to the transferred amount of charge, is gen¬ 
erated. This internal field breaks the inversion symmetry 
and thus generates RSOC in the respective layer. 

It stands to reason that also the electrons in the sur¬ 
face layer LI experience a RSOC because the same elec¬ 
tric field acts on the electronic states in LI. In the case 
of a bare RSOC 7 = /XR/SmeC in 7 (E x p) • S for free 
electrons with mass We, momentum p and spin S in a 
plane subject to a transverse electrical field E, a reason¬ 
ing demanding an equivalent treatment of layers LO and 
LI with respect to the RSOC appears to be appropri¬ 
ate. However, the bare RSOC with electric fields gen¬ 
erated by the charged layers is orders of magnitude too 
small to have sizeable impact on the electronic recon¬ 
struction in the heterostructures. The resolution of this 
apparent unfeasibility of bare RSOC effects has already 
been established for semiconductor heterostructures (see 
the book by R. Winkler [15] for a comprehensive treat¬ 
ment) and is implemented for the LAO/STO interfaces 
as follows m IS]: The Ti-ions provide an atomic SOC, 
whereas an internal electric field which displaces the 
Ti-ions slightly out of their in-plane positions, couples 
the three Og orbitals. As spin and orbital space cou¬ 
ple, one has to diagonalize six bands. The diagonaliza- 
tion of the 6 -band system generates a SOC of Rashba 
type ittRfc (see Eq. 0 ) even though only the lowest 
of the Og bands (the dxy) may be partially occupied. 


The RSOC parameter or has to be estimated from elec¬ 
tronic structure evaluations [5]. In fact, from magne¬ 
totransport measurements, a sizeable RSOC has been 
identified for the LAO/STO interface electronic state 
(ag= « (1...5) X 10-i°eVcm « (0.74... 3.7) x IQ-^e^ 
for different backgate voltages), where the value of 
depends on the charge carrier density which is controlled 
by the internal electric field. In Sec. jlVj we will address 
the 6 -band situation explicitly. Until then we consider 
an effective two-band situation with a RSOC that is con¬ 
trolled by the internal field that depends in turn on the 
charge density. The actual set-up for the electronic sur¬ 
face states (layer LI) is very different and SOC effects 
have so far not been found to be relevant. Therefore we 
disregard RSOC in the surface layer LI in our approach. 

Through the mechanism outlined above, the strength 
of the induced RSOC, or, becomes density dependent: 
We denote the Helmholtz free energies of the single layers 
as Fo{n) and Fi(ni), calculated in the first stage without 
connection between LO and LI. In fact, we calculate 
the single layer grand canonical energies and then use 
a Legendre transformation to obtain the Helmholtz free 
energies. The total free energy of the connected system 

Ttot = Fo{n) + Fi{ni) + Fefin^ni) 


is the sum of the single layer energies and the electro¬ 
static term Egg. The total energy Ft^t is a function of 
the total density Q = n + ni only, that is, the sum of 
both single layer densities: the dependence on the single 
layer densities is eliminated by the condition that n has 
to minimize the total free energy, which is equivalent to 
the equality of the electrochemical potentials of both lay¬ 
ers. Note that, as the layers are connected, or becomes 
a function of the electric field generated by LI and thus 
of Q and n{Q). 

Since both layers are coupled through the electric field 
and are supposed to be in electrical contact, only the 
inverse compressibility of the total system 


K 


-1 




d2 (F/A) 

dg2 


( 2 ) 


is a well-defined thermodynamic quantity. Here, 
F = F{Q,V) is an appropriate Legendre transform of 
Ftot{Q,n) with respect to the interface charge density 
n: the thermodynamic variable is then the voltage bias 
V between interface and surface (see Appendix [d| . The 
electronic compressibility derives from the second deriva¬ 
tive of the electronic energy with respect to the volume. 
In the case of two-dimensional electron systems (as in the 
considered heterostructure) it is the area instead of the 
volume. This translates into the second derivative with 
respect to the electronic areal density Q. The second 
derivative in Eq. ([^ is to be taken with respect to the 
total dependence on Q, that is, also the Q dependence of 
the internal field in the polar film enters. The positivity 
of K then represents a stability criterion with respect to 
charge separation in the plane directions. 
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FIG. 2. (Color online) (a) Multilayered system with static (induced) charges (red/blue), electrical potential perpendicular to 
the layers (orange) and electric fields (green), (b) Electrostatic equivalent system which relates to the configurational build-up 
of a polar LAO-film on top of an LAO/STO interface layer. 


III. POLAR MATERIAL WITH A RSOC 
INTERFACE TO A NON-POLAR INSULATOR 

To approach a polar system with a surface and an 
interface to a non-polar insulator, we first present a 
model system that consists of three two-dimensional lay¬ 
ers (Fig. I^a)) but only two of them comprise mobile 
electrons: The planes at the bottom and at the top have 
areal charge densities of —ne and —{Q — n)e, respec¬ 
tively. The electrons may be transferred between these 
two planes to minimize the free energy, while the plane 
sandwiched between them—with a distance di to the bot¬ 
tom and d 2 to the top—holds immobile, fixed charges of 
density Q. This charge density Q compensates the charge 
of the other two planes and therefore keeps the system 
neutral. The dispersion of the electrons in the bottom 
layer (LO) is given by Eq. ([^, while the top layer (LI) is 
described as an electron system with a density of states 
(DOS) of iVi. If n ^ 0, an electric field perpendicular to 
LO exists which derives from the total charge of the two 
other layers. This field induces the RSOC in the bottom 
layer and is proportional to the electron density n of LO. 
The electrostatic energy per area of this layout is 

u. 27rdie2^2 , 27rd2e%^ ^^2 

Fes(0,u)/A — 7Z -t“ (0 u) , (3) 

£1 £2 

where are the dielectric constants of the material 
between the planes (see Appendix E- This arrange¬ 
ment is an effective model for the oxide heterostructure 
LAO/STO, where polar layers with a distance d{ and 
^2 between them alternate (Fig. |^b)). Due to elec¬ 
tronic reconstruction, electrons from the surface relocate 
to the former empty layer at the LAO/STO interface 
and occupy Tis^^ states [20l|2T]. We assume that the ef¬ 
fective distances to the sandwiched layer are equal, i.e., 
di/ei = d 2/£2 = d/{2e). In LAO, the distances di and 


d 2 are approximately equal but £i /£2 is not accurately 
known. Therefore we take a layer-independent dielectric 
constant e for the LAO film. This simplification is of no 
qualitative consequences for the further evaluation; the 
results depend on an effective distance d/e. 

There are various models for the density dependence of 
au (cf. Ref. [H]). In general, the strength of the RSOC 
au is taken to be proportional to the SBEF: or oc E. 
The SBEE E is induced by the charge carrier density -l-ne 
(see Fig.J^a)) and can be expanded as E{n) = Eq + £ ■ 
n + 0{n^, where £ ■ n = Ei/2. The factor 1/2 accounts 
for the fact that the SBEF is generated by the charge 
carrier density +ne at the sandwiched layer and not by 
the electron density —ne at the LAO/STO interface LO. 
The field Eq refers to an electric field that is already 
present in the absence of a charge density +ne\ it may be 
produced by dipolar distortions in the LAO-layer next to 
LO. With this expansion of E{n), one can also decompose 

aR(n) = Qfo + cki • n-b O(n^). (4) 

A more adequate determination of Q;R(n) will be given 
in the context of the 6-band model for the LAO/STO 
interface in Sec. EYl 

In order to implement the difference in chemical po¬ 
tentials of the electronic systems at the interface LO and 
surface LI, we introduce a term —Vqu in the total free 
energy Ftot: 

Ftot{Q,n) = Fo(n) + Fi{Q - n) + Fes{Q,n) - cVq nA. 

(5) 

The potential Vq may be determined from electronic 
structure evaluations; we use it phenomenologically to 
tune the electron density at layer LO to the experimen¬ 
tally observed values. 

The thermodynamic relation dF^ot/dn = eA V{Q, n) is 
used to find the function n{Q, V) (here we still consider 
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FIG. 3. (Color online) (a) Inverse compressibility over plate distance dje and strength of RSOC a = ao + ai n. The blue plane 
indicates = 0. Here we take mo = 0.7me, mi = me, Q = le“/uc = 6.6 • 10^^/cm^, Vo = 0, and ao = 0. (b) Phase of the 
system: PI"*" (yellow), P2 (green), Plo (cyan) and P2o (blue). 


V = 0, and n{Q,0) = n{Q)). Since n{Q) is limited to 
the interval [0, Q\, it is possible that the boundary values 
n = Q or n = 0 yield the global minimum oi F{Q,V). If 
neither all nor no charge resides at the bottom layer, an 
analytical expression for the stability of the total system 
can be derived (Appendix [d|) : 


K > 0 


dn 


< 1 . 


( 6 ) 


We would like to point out that a system with /c < 0 
can still be stabilized by a mechanism not related to the 
electronic system in LO and LI (and thereby yield a pos¬ 
itive total compressibility). In this case, an enhancement 
of the total electronic charge in LI causes an amount of 
charge larger than the added charge to flow to the RSOC 
layer LO. The distribution of charge between LO and LI 
is controlled by the electric field which is in turn deter¬ 
mined by the positive capacitance of the system. On the 
contrary, if the system cannot be stabilized (that is, the 
total compressibility is negative) then it is expected to 
develop a phase separated (inhomogeneous) state in the 
plane. This latter scenario will not be investigated here 
but we will comment in our conclusions on the proposi¬ 
tion EHII] of RSOC-driven phase separation. 

We will not discuss the temperature dependence of the 
capacitance and therefore reduce our considerations to 
zero temperature in the following evaluations. The free 
energy per area of the top two-dimensional layer (LI) 


Fo{Q-n)/A=^JQ-nf (7) 

depends on the electron density transferred to the bot¬ 
tom layer. Here we assumed an interactionless electron 
system in LI, characterized by a density of states Ni. 
With Ni —>■ oo (that is, the mass of the charge carri¬ 
ers in LI increases beyond limits, mi —>■ oo), a localized 


state can be generated technically. The surface of the 
LAO/STO heterostructure is insulating. Then, an addi¬ 
tional metallic layer on top of the heterostructure pro¬ 
duces a capacitor. The top electrode can be a gold or 
YBa 2 Cu 307 (YBCO) film, which connects well with the 
topmost AIO 3 layer (see Appendix]^. There, the capac¬ 
itance increase is slightly smaller than for a 2D electrode 
due to the finite screening length in the 3D film but the 
qualitative results are the same as presented in this sec¬ 
tion for the 2D top electrode. 

We re-emphasize that the top layer LI must be in¬ 
cluded necessarily in order to render the correct electro¬ 
statics which then yields the proper total compressibil¬ 
ity and capacitance. However, as long as interaction or 
spin-orbit effects are negligible in LI, the impact of LI 
is determined by two parameters, only: the density of 
states Ni at the surface electrode and the energetic dis¬ 
placement —eVb between interface and surface electronic 
levels. The latter controls trivially the electronic density 
n at the interface. The former yields the kinetic term 
in the quantum capacitance and is known to produce 
a slightly increased effective distance between the elec¬ 
trodes (see, e.g., Hefs. IT^ llBl [22U24)l . Consequently, the 
materials properties of the surface electrode do not af¬ 
fect the discussed physical nature of the heterostructure 
qualitatively, except for adjusting the carrier density at 
the interface. 

The energy per area of the layer with RSOC (with 
apt as density independent constant), calculated in Ap¬ 
pendix |AT] is 


F^{n)/A 


H? TT 2 

ri^ 2 ,r3 4 

—- 'KNoa^n -\ -Yifat 


< 0 
,^l>0 

( 8 ) 


When taking the density dependence of api into account. 
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we expand aR(n) in n (cf. Eq. Q) and keep terms up to 
linear order in n: 

a = ao + ai ■ n. (9) 

One has to distinguish between the case where only the 
lower band of the RSOC split-bands is occupied, (see 
Fig. [^with n < Ob and the case where both bands are oc¬ 
cupied, (Fig.mwith > 0). Since ^ = 0 corresponds 
to n = TrfVga^, the first case equals either n < n_ or 
n > n_|_, with 

= 2 ^(jVoai)2 ■ 2^7V2aoai ± - dTriV^aoOi) . 

( 10 ) 

This means that also for large densities n > n+, all 
electrons are in the lower band, because the minimum 
energy of the bands, = —7J‘.^o(ao + ain)^/2, 

decreases for increasing density. 


one electron per areal unit cell of the LAO/STO struc¬ 
ture. The transition to a regime of negative compressibil¬ 
ity is indicated by the transparent blue plane: the area 
below the blue plane exhibits negative compressibility, 
which can occur for all phases except PO. For phase P2, 
however, the negative compressibility regime is small, re¬ 
stricted to very small d/e (less than ~ 0.2 ob) and large 
RSOC («! > lO(eaB)^). For large effective plate dis¬ 
tances both layers are always occupied, since the electro¬ 
static term is then dominating. In the reverse case—for 
small very d/e —the electrostatic cost for transferring all 
electrons to one layer is negligible. Then, as n+ decreases 
with increasing oi, only one band in LO is occupied for 
strong RSOC. 

For Oo = 0 the phase boundary between P1+ and PIq 
is given by 

ire^d {SttNqoIQ)^ — 1 
^ UnN^alQ ' 


With the total energy given by Eqs. and (|^ , we 

calculate the minimizing density n{Q) numerically, and 
analytically for oq = 0. We identify six different solutions 
for n{Q) and determine the respective compressibility: 
n = 0 (PO): The bottom layer LO is empty and the total 
compressibility is given by the compressibility of the top 
layer. 

0 < n < n_ (Pl“): This phase exists only for ao > 0. 
All electrons of the bottom layer are in the lower band. 
The compressibility is calculated numerically according 
to Eq. (iD^. 


n_ < n < min((5,n_|_) (P2): Both bands of the bottom 
layer are occupied. The co mpres sibility is calculated nu¬ 
merically according to Eq. (D5b). 


= Q < n+ (P2o): All electrons reside in the bottom 


layer, and both bands are occupied. If n = Q is inserted 
into the total energy Ftot and the second derivative with 
respect to Q is taken, the inverse compressibility is 


/ 

«P2o ^ “ 27rAroQ;i(2ao -l- SaiQ) 

+ 27r^Afgai(ao -|- aiQ)^ -I- 27re^A^o^ 


( 11 ) 


n+ < n < Q (PI'*'): All electrons of the bottom layer 
are in the lower band. The compressibility is calculated 
numerically according to Eq. (D5b I. 
n+ < n = Q (Plo): All electrons reside in the bottom 
layer, but only the lower band is occupied. Inserting 
n = Q into Ftot yields for the total compressibility: 


-1 f QqQ 

No \TTN^{ao + aiQ)^ 

— ttA^q Oi (2 q!o- l-SoiQ)-|-27re^ A^g —^ . (12) 

In Fig. I^a) we display the inverse compressibility in the 
respective phases as a function of the effective plate dis¬ 
tance d/e and of ai for oq = 0 and a total charge of 


This relation can be used to show analytically that kpiq 
is always negative. One deduces from Eq. (0 for phase 
P2o that also kp 2^ < 0. Obviously, the phases PIq 
and P2o display a “phase separation between the lay¬ 
ers” (the index 0 indicates that all mobile charge carriers 
are in layer LO). Moreover, the negative compressibil¬ 
ities (kpio < 0 , Kp 2 g < 0) signal that these phases are 
thermodynamically unstable and that a phase separation 
within plane LO is expected. 

When no charge resides on the surface—as it is the 
case in phases Pig and P2o— one could expected that 
the compressibilities of Eqs. © or ( [I^ are identical to 
the compressibilities calculated for one isolated layer in 
Ref. [9] (apart from the contribution of the electrostatic 
energy which straightforwardly produces the geometrical 
inverse capacitance). This is, however, not the case: our 
results differ from those of Ref. [9], as explained in Ap¬ 
pendix 1^ This discrepancy can be traced back to the 
inclusion of a n-dependent RSOC aR(n) either directly 
in the free energy or in a later stage in the chemical po¬ 
tential. We calculate the inverse compressibility through 
the second derivative of the free energy with respect to 
Q = n whereas the authors of Ref. [S] identify the inverse 
compressibility through a first derivative of the chemical 
potential. The latter approach misses a non negligible 
contribution from the derivative of the energy with re¬ 
spect to the density in aR(n). For the phase Pig, for 
example, the latter approach results in an inverse com¬ 
pressibility that is three times smaller than our result for 
phase Pig. Moreover, those phases which do not have all 
electrons transferred from the surface LI to the interface 
LO, cannot be addressed by the single layer scheme in 
Ref. [9] , and we identify them from the free energy which 
comprises an n-dependent RSOC aR (n). 

The application of a voltage bias V between interface 
and surface, i.e., between LO and LI, can be used to tune 
the electron density on the interface and hence also the 
phase of the system (see Fig. [^: As n increases, so does 
the strength of the RSOC a, and the band minimum is 
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FIG. 4. (Color online) Inverse compressibility, areal density, and capacitance for a two-band interface model as a function of 
the bias voltage V, for which e/as corresponds to 27.2 V in SI units. Here mo = 0.7me, mi = me, d/e = IOab, Q ~ le“/uc, 
ao = 0, and ai = 20e^aB for a two layer system, (a) Inverse compressibility of the total system obtained as analytical solution 
for the one-band case (blue). Density n on the interface layer (red), (b) Differential capacitance per area calculated analytically 
(blue), geometric capacitance per area (black dashed), capacitance per area without RSOC (brown dashed) and capacitance 
eAn/Ay per area (orange). The background color indicates the phase of the system according to the caption of Fig. while 
violet corresponds to PO. 


shifted to lower energies until only one band is partially 
occupied and the system is in state P1+. This process 
increases the compressibility that stays positive even if 
all charge is transfered to the interface (Plo)- This does 
not contradict the analytical analysis from above, which 
was for V = 0. 

As for the phases PO and PIq, note that kpq is inde¬ 
pendent of ao, a I and Nq, as no charge is in LO, and that 
kpiq does not depend on the second electrostatic term in 
Eq. (§, as no charge is in LI. This observation explains 
the sharp transitions in k{V) in Fig. |^a). 

The differential capacitance Cdiff is derived in ap¬ 
pendix 1^ 


A/Cdiff = 


1 W 

e dn 


dlFp + + ggfes 

A ■ 


(13) 


The term d^F^g yields the geometrical capacitance 
Cgeo = Ae/Aird. The total capacitance is enhanced 
above the geometrical value, if the sum of the inverse 
capacitances of the single layers is negative. In Fig. 
we compare Cdis, the capacitance C = e^njtFV^ the 
geometric capacitance Cgeo, and Cq=o = A/{^'Kd/e + 
1 /e^Ao + 1/e^A^i), the capacitance for an interface with¬ 
out RSOC. It is evident that the differential capacitance 
is enhanced above the capacitance of a system without 
RSOC for finite interface electron densities. Moreover, 
both Cdiff and C increase beyond Cgeo for sufficiently 
large V (Fig.|4];b)). However, this does not signify a neg¬ 
ative compressibility k (cf.Fig. |^a)) —although the par¬ 
tial compressibility of the interface my be negative (the 


latter is not a well-defined thermodynamic quantity). As 
elaborated in Appendix the connection between the 
compressibility of the system and its differential capaci¬ 
tance is 


A = Q' 


.d^F 

W 


( Cdiff 

e2A2 aQ2 ) ’ 


(14) 


where F{Q,V) is the Legendre transform of Ftot{Q,n): 
F{Q, V) = FtotiQ, n{Q, V)) - eA P • n(Q, P). With the 
general assumptions that the total energy Ftot of the sys¬ 
tem is a functional of n and Q — n exclusively, and that n 
is not a boundary value {Q or 0) then the compressibil¬ 
ity of the system is related to the differential capacitance 
through 


1 e^A / dn\ dn 


(15) 


where both, k and Cdiff, are functions of Q and P. The 
first factor on the rhs is what one should expect for a 
single metallic electrode: the electronic compressibility 
is the capacitance of the electrode up to a factor e^Q'^A. 
The further two factors with the derivative dn/dQ ac¬ 
count for the charge redistribution between surface and 
interface—which, of course, depends implicitly on the 
RSOC (see Appendix [pj). 

Finally, we discuss the size of the experimental values 
of the control parameters. We use mo = 0.7me for the 
effective mass of the electrons at the LAO/STO interface. 
This value results from the multi-band evaluation [1] , ad¬ 
dressed in the next section: a light mass m^ = 0.7 nie and 






















FIG. 5. (Color online) Energy bands for m{ = 0.7me, = 15me, Aso = lOmeV, Ae = 50meV, Az,o = 0.61 • 10 ^ e^, and 
Az,i = 8.29 e^ag. The bands are plotted from the F-Point to the A-point (continuous curves) and to the M-point (dashed 
curves). The highest occupied energy level for the densities n = 2 ■ 10^^/err? (a) and n = 2 • 10^^/cm^ (b) is indicated by the 
brown curve. 


a heavy mass mu = 15 rrie were identified from ARPES 
measurements [25j for the electrons of the dxy-, dxz- and 
dyz-bands. Joshua et al. illustrated in the Suppl. of 
Ref. |1] that, within the multiband model, the electrons 
of the lowest resulting band acquire an effective mass of 
Too = 0.7TOe close to the band minimum. 

Estimating ao and ai has the largest uncertainty: 
Caviglia et al. [3] concluded from their magnetotrans¬ 
port experiments that ag' « (1...5) x 10“^^eVm « 
(0.74... 3.7) X 10“^ in the range of backgate volt¬ 
ages where the superconducting dome is formed. For 
a more negative bias the interface becomes depleted 
and the RSOC is nearly constant with a low value of 
ao ~ 10“^^eVm « 0.74 x 10“^ e^. We neglect this 
rather small constant contribution with respect to the 
term linear in the bias. The RSOC increases approxi¬ 
mately linearly with bias in between values correspond¬ 
ing to the low-density foot and top of the superconduct¬ 
ing dome (compare Fig. 3 in Ref. [3]). The difference in 
charge carrier density at foot and top of the dome [55] is 
« 0.015 e“/uc = 9.9 x 10^^/cm^, so that 


4 • 10“^^ eVm 
9.9 • 10^^ cm“^ 
lO(eaB)^. 


4.0- 10"^^ eV-cm^ 


(16) 


We do not know the exact values for the density at which 
aff was measured. Therefore ai « 10 (cob)^, as deduced 
from the experiment of Caviglia et al. |3] , might be some¬ 
what smaller but certainly not more than one order of 
magnitude. 


IV. MULTIBAND LAYER WITH SPIN-ORBIT 
COUPLING 


The free electrons at the LAO/STO interface reside 
predominantly in the titanium 3d orbitals of the first 
TiO layer. For a refined description of the interface 
we use a six-band model for the three spin-split t 2 g 
orbitals m 0 Ez). Within this formalism, the emer¬ 
gence of a strong Rashba-like SOC is explained naturally 
through the breaking of inversion symmetry at the inter¬ 
face, which couples the atomic SOC of the Ti atoms to 
the inter-band hopping. The Hamiltonian can be written 


in the dy^, d^x, dxy basis as: 


H = Hk + Hso + Td-z. 


(17) 


The first term is the intra-orbital hopping: 

/ 




277 A^h 


k^ 

277 A^i 


V 


0 

I feg 

277iV’i ' 277A^h 

0 


0 
0 


>CTo, 


(18) 


Ni = mi/{Trf?) is the density of states for light (mi) or 
heavy (wh) electrons and Ae the splitting between the 
dxy- and the dy^-, dzx-orbitals. The contribution from 
atomic spin-orbit coupling is 

diso = Aso • ^ Li C) cTj 


(19) 
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/h'l 0 0.5 1.0 1.5 2.0 

V{e/a,) 


FIG. 6. (Color online) Inverse compressibility, areal density, and capacitance for a six-band interface model as a function of 
the bias voltage V. Here mi = 0.7 me, m^ — ISrrie, mi = me, dje = IOub, Q = le“/uc, and Aso = lOmeV, and e/oB 
corresponds to 27.2 V in SI units, (a) Inverse compressibility of the total system (blue) and density on the interface (red) for 
Az.o = 0.61 • 10“^ e^, Az,i = 8.29e^a|, As = SOmeV, and Vb = —1.149e/aB. The background colors indicate the number 
of occupied bands: two (yellow), four (green) and five (purple). The brown circles refer to systems with bias values V that 
correspond to the densities n of Figs. a) and (b), respectively, (b) Differential capacitance per area for three different sets 
of parameters with n = 2 • 10^^/cm^ at H = 0, which is provided by Vb = —1.149e/an (blue), Vb = —1.149e/oB (green) 
and Vb = —1.155e/as (red). For small k, the energy splitting Ai5 (as in Fig. for the blue and red curve corresponds 
to ai = 2.5e^a| and for the green curve to a\ = 5.0e^a|. For the dotted part of the curves, a fc-space region around the 
k = 0 wave vector is unoccupied. The geometric capacitance per area (black) and the differential capacitance per area for 
Vb = —0.781 e/aB, Az_o = 0.3 • 10“® and A^^i = 0 (brown) are also shown. 


/ /o 0 0\ /o 0 -A / 0 i 0\ \ 

L= 0 0 f , 00 OL -fOo bi (20) 

\\0 -i Oj \i 0 0 J \0 0 oj J 

where Aso is the strength of the atomic SOC, ao is the 
2x2 identity matrix, and ai are the 2x2 Pauli-matrices. 
Due to the confining potential at the interface and its 
asymmetry, the d-orbitals are deformed and hence no 
longer orthogonal. This allows (anti-symmetric) interor¬ 
bital hopping between dxy- and dxz-orbitals along the y- 
direction and between dxy- and dyz-orbitals along the x- 
direction, presented by 

/ 0 0 —*fcx\ 

•H, = Aj 0 0 iky \ OcTo- (21) 

yifcx —iky 0 / 

We assume that the magnitude of the hopping is pro¬ 
portional to the asymmetry inducing electrical field. In 
our modeling, the origin of this electric field is the (posi¬ 
tive) charge on the surface that compensates for the elec¬ 
trons transferred from the surface to the interface system. 
Therefore we write: 

Az = Az 0 + Az 1 • n. ( 22 ) 

The term proportional to n is generated through 
the asymmetry inducing electrical field, and the n- 


independent term Az_o can be traced to dipolar distor¬ 
tions in the LAO-layer next to the interface. The di- 
agonalization, which can be done analytically only for 
kx = ky = 0 , yields three spin split doublets. Hz, which 
originates from the deformation, vanishes at k = 0 and 
the doubly degenerate energies are 

ei(fcx = 0, fcy = 0) = ^ (-Ae + Aso - S) , (23a) 

e 2 (fcx = 0 , fcy = 0 ) = -Aso, (23b) 

£3{kx — ky = 0) = - (—Ae + Aso + S) (23c) 

with S = \J A| -I- 2AeAso + 9^io- 

For a non-zero Az, the degeneracies are lifted. Moreover, 
crossings between different doublets are avoided [1 [g 
(see Fig. [^). With increasing A^, the splitting between 
the formerly degenerated bands becomes larger and the 
point of lowest energy in fc-space moves to larger values 
of k (cf. Fig. [^b)). 

The shape of the bands depends, analogously to the 
two-band model with RSOC (cf. Eq. Q), on the electron 
density. We add a top gate as in Sec. |III| and perform the 
same minimization of the energy of the combined system 
with respect to n for different applied voltages. The re¬ 
sults are shown in Fig. where we chose the light and 
heavy electron masses = 0.7me and mh = 15me, 
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respectively (in accordance with the ARPES mea¬ 
surements for SrTiOa surfaces). The energy splitting is 
taken to be Ag = 50meV, according to the x-ray mea¬ 
surements in Ref. [55]. Alternatively, we also consider 
Ae = 150 meV, which is in better agreement with vari¬ 
ous ab initio evaluations (see, e.g.. Ref. 0 ). The atomic 
SOC was estimated from transport experiments 11110 
and ab-initio calculations m to be Aso = 10-25 meV; 
we chose the lower boundary of 10 meV. 

The ^-dependent energy splitting AE{k) induced by 
Az in the lower and upper doublet increases linearly 
for small k, i.e., it is Rashba-like: AE{k) = a(Az)fc 
(cf. Fig.0. This can be employed to fit the parameters 
Az^o and Az_i to the estimates of oq and ai made from 
experimental results in the previous chapter: by writing 

a{n) = a{A^{nj) = 6 (Az,o + A^j • n), (24) 

we identify bA^^o = «o and bA^^i = ai, where b is 
obtained from fitting a(Az) to the numerically calcu¬ 
lated AE{k). For Ae = 50meV and Aso = lOmeV 
this yields b « 3.3 and we find A^ g ~ 2.5 • 10“^e^ and 
Az^i « 33e^a|. Although this estimation for A^^g and 
Az^i neglects the dependence of AE on the angular di¬ 
rection in /c-space, it gives a reasonable approximation 
for a{n) as long as the charge density is small, such that 
AE{k) oc k holds. For larger densities, for which this pro¬ 
portionality fails around the Fermi momentum (i.e., it is 
not Rashba-like, cf. Fig. [^, the above estimation is no 



k ) 

FIG. 7. (Color online) Energy splitting between the two 
lowest bands (blue), the third and fourth band (green) and 
the two highest bands (red). Here the parameters are 
rai = 0.7rue, rrih = 15me, Aso = lOmeV, As = 50meV, 
Az,o = 0.61-10“^ e^ Az,i = 8.29 e^a|, and n = 2 • 10^'‘/cm^ 
as for the right side of Fig. The energy splitting is plotted 
from the F-Point to the A-point (continuous curves) and to 
the M-point (dashed curves). 


longer valid. The capacitance and compressibility of the 
system is given mostly by the electrons from the lowest 
doublet, because for low densities the higher bands are 
unoccupied or considerably less occupied than the lowest 
two bands. For larger densities, on the other hand, the 
variation of the electron number upon changing V or Q 
is largest in the lowest doublet (since then the DOS is 
larger in the lower doublet, cf. Fig. [^b)). 

The differential capacitance Cdis for different gate 
voltages V is displayed in Fig. |^b). The difference of 
the chemical potentials in the interface and the metallic 
top layer, Vg, was set to a value that gives a charge den¬ 
sity n = 2 X 10^^/cm^ for V = 0. For A^^i = 0, i.e., a 
density independent bandstructure, the capacitance in¬ 
creases with voltage (see brown line in Fig. [^b)) until 
the Fermi energy reaches a value where the deviation 
from a parabolic band structure becomes negligible and, 
consequently, the DOS is constant. Notice that the ca¬ 
pacitance derived from the microscopic model is smaller 
than the geometric capacitance, since the kinetic energy 
of the charge carriers in the electrodes adds a positive 
term to the inverse capacitance (see, e.g.. Ref. m)- 

A non-zero A^ i enhances Cdis, even above the geo¬ 
metrical capacitance if A^ i is sufficiently large. With 
increasing voltage V, the density n at the interface in¬ 
creases (Fig. [^a)), and hence A^. This flattens the dis¬ 
persion of the lower doublet and thereby enhances the ef¬ 
fective mass at kp (see Fig. [^, which causes the increase 
of capacitance above the value for A^^ = 0 (cf. Figl^b)). 
For large A^, the energy minima of the two lowest bands 
are far from the zero point and k = 0 is not occupied. In 
this regime, AE{k) is no longer linear in k and the esti¬ 
mations for Az_g and A^ i are no longer valid. This region 
is indicated by the dots of the differential capacitance in 
Fig.pb). 

Within the regime where the bands become flatter with 
increasing density, negative compressibility may appear. 
The inverse compressibility for different values of A^ i is 
shown in Fig. 0 as a function of distance between the in¬ 
terface and surface. Evidently, Fig. [^a), which presents 
the inverse compressibility for the 6 -band model corre¬ 
sponds to Fig. I^a) for the 2-band model. For a small 
effective distance d/e, all charge resides on the interface 
(Fig. I^b)), the compressibility is negative, and k~^ in¬ 
creases nearly linearly with the distance. This is under¬ 
stood by setting n = Q in Eq. (|^ —the only term in 
the total energy with a dependence on d/e —and taking 
the total second derivative with respect to Q. For large 
distance the electrostatic energy dominates and the to¬ 
tal energy is minimized by n « Q/2. Inserting this into 
Eq. ® yields half the gradient of the n = Q regime, see 
Fig. ^a). 

In Fig. I^a) the inverse compressibility is displayed for 
Vg = 0, whereas in Fig.|^c) the potential is Vg = —31 eV. 
For Vg = 0 and Aso = 0, the lower band edge at the in¬ 
terface is shifted by — Ae = —50meV with respect to 
the band edge of the surface electrode. Then, the elec¬ 
trons tend to go to the interface if no other effect such 
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FIG. 8. (Color online) Inverse compressibility for Vo = 0 (a) and for Vo = —31V (c), both as functions of the effective 
interface-surface distance d/e. Here mt = 0.7rrie, rrih = 15me, mi = rrie, Q = le“/uc, Aso = lOmeV, Ae = SOmeV, and 
Az,o = 10“^ e^. The plots below, i.e., panels (b) and (d) display the interface electronic density n for the respective potentials 
Ho.’ 


as the RSOC interferes, provided that the electrostatic 
energy is small—and the compressibility can be negative 
for sufficiently large A^^i. 

However, the behavior of as a function of d/e 
changes drastically when we adjust Vb to —31V so that 
n = 2 X 10^^/cm^ for d/e = 10 ob and = 33e^a|. 
Then n is in the range of carrier densities established for 
the LAO/STO interface with d/e ~ 10ob- The inverse 
compressibility first increases linearly with d/e because 
n = 0 in the regime of small distances (see Fig. ic)). 
There is a small offset due to the finite compressibility 
of the surface electronic state. Then, with a further rise 
in d/e, the polar-catastrophe mechanism induces a fi¬ 
nite electron density at the interface (Fig. I^d)), and the 
inverse compressibility decreases sharply ^ig. [^c)) be¬ 
cause a small but finite n reduces the electrostatic term 
linearly with n (see Eq. §)■ Eventually, k ^ increases 
towards its limiting linear behavior for large d/e as in 
Fig. [^a). This analysis allows us to conclude that the 
compressibility stays strictly positive even for large val¬ 
ues of the RSOC if the system parameters are adjusted 
within a range that is consistent with parameter values 
previously identified for the LAO/STO heterostructure. 

For a thorough determination of the compressibility of 
the LAO/STO heterostructure additional effects may be 
relevant: these include oxygen vacancies, their distribu¬ 
tion between the planes as well as the lateral extension 
of the two-dimensional electron liquid at the interface. 


However, these effects influence primarily the charge car¬ 
rier density. As the charge carrier density is adjusted to 
the experimental values through a bias Vq in our analy¬ 
sis, we expect that the evaluation of the compressibility 
remains valid. 


V. CONCLUSIONS 

Low-dimensional charge states at interfaces of complex 
oxides have been explored intensively in recent years, and 
it is a combination of scientific and technological inter¬ 
est which drives the research in this field. This is also 
valid for the work presented in this article: we explored 
the electronic charge redistribution in a heterostructure 
of a polar and a non-polar insulator, the paradigmatic 
realization of which has become the LAO/STO interface 
of a polar LAO film on the non-polar substrate STO. 
Understanding the quantum phenomena that control the 
electronic reconstruction at the interface, and thereby 
also the charge redistribution between surface and in¬ 
terface, may provide a possible strategy to control ca¬ 
pacitances by the electronic properties of the conduct¬ 
ing electrodes—as opposed to the standard search for di¬ 
electrics with high dielectric constant [3D]. 

Sizeable capacitance enhancements and negative 
electronic compressibility have been observed for 
LAO/STO [HID]. Recently, these findings were related 
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to a strong RSOC [5HTT]. The RSOC is generated by 
an electric field that in turn is controlled by the elec¬ 
tronic charge density —en at the interface. In fact, a 
band structure, which depends on the density of the band 
electrons themselves, allows for a decreasing chemical po¬ 
tential with increasing charge carrier density—assuming 
that the RSOC strength is sufficiently large. Such a de¬ 
pendence of the chemical potential on the carrier den¬ 
sity represents a negative compressibility. Thus, instead 
of the standard explanation for a negative compressibil¬ 
ity in low density semiconductor systems m or in oxide 
heterostructures mis] through exchange and correlation 
effects, it may be that spin-orbit effects play an essential 
role. 

However, in order to settle the realization of a negative 
compressibility and the possible formation of an inhomo¬ 
geneous charge-separated state, one has to investigate at 
least the compressibility of the complete electronic sys¬ 
tem (in the case that all electronic subsystems are cou¬ 
pled) . The compressibility of other subsystems which are 
not comprised of mobile electrons may also contribute 
and keep the total compressibility positive. In this arti¬ 
cle we not only included the mobile interface electrons, 
but also the surface charges with their respective field. 
This is to be considered as a minimal setup for two rea¬ 
sons: (1) the RSOC is controlled by the electric field of 
the charges at the surface, and (2) this electric field and 
the charge carriers of the surface electrode contribute to 
the compressibility of the interface-surface system. 

A remarkable relation, which only involves the capaci¬ 
tance of the interface itself, Cq, and the geometric capac¬ 
itance Cgeom = eA/And of the interface-surface system, 
allows to determine if this system attains a negative com¬ 
pressibility: 

K < 0 ^ + 2^ <0, (25) 

-A S 


where Cq^ = (A/eA)"^ is found from the sec¬ 

ond derivative of the interface free energy Fq with re¬ 
spect to the interface electronic density n. This relation 
was derived in Appendix [Plunder the conditions that the 
interface-surface capacitance Cjjg = (I/eA)^ d'^F/din? is 
positive and that n does not assume a boundary value 
(that is, neither n = 0 nor n — Q holds). Here, F is 
the free energy of the entire heterostructure. The rela¬ 
tion (251 is remarkable in two respects: first, the com¬ 


pressibility K can become negative even when the capac¬ 
itance Cdifif is positive, and second, this stability crite¬ 
rion does not depend explicitly on the compressibility of 
the surface—even though the total inverse compressibil¬ 
ity K~ ^ contains a contribution from the surface electron 
system. Both statements can be traced back to the spe¬ 
cial interface-surface setup in which the capacitance is 
not equal to the total compressibility (see Eq. (HI)- 
The relation (|25|) can be cast into a more intuitive form 


(Appendix : 


n 1 


This inequality obviously expresses the fact that the sys¬ 
tem is stable if an additional overall electronic charge 
6Q does not induce a charge transfer Sn larger than 
the added overall charge. The compressibility is nega¬ 
tive and the electronic state may charge-separate only if 
dnjdQ > I. 

In order to be explicit, we considered first a 2-b and 
model for the interface electronic system (Secs. HI and 
Appendix |C|). The bands are split by a RSOC that de¬ 
pends on the charge density of the surface electronic 
system. The charge densities at interface and surface 
were calculated from the minimization of the total free 
energy. This simple model allows to identify a series 
of electronic phases that depend on the occupation of 
the interface split bands. The second derivatives of the 
free energy with respect to the total electronic density 
Q and the interface electronic density n yield (phase- 
dependent) compressibility and capacitance, respectively. 
A negative compressibility arises only for very small ef¬ 
fective distances d/e between interface and surface, and 
strong RSOC (see Fig. I^a); for example for an RSOC- 
parameter ai ^ 10 (euBj the effective distance has to be 
d/e < 0.2Ob ). Even for RSOC parameters twice as large 
as the largest RSOC estimated from magneto transport 
measurements for LAO/STO, the compressibility stays 
positive (cf. Fig. |^a)), also for large but sensible val¬ 
ues of a bias voltage V. This differs from the findings 
of Refs. pH [TH]. The discrepancy originates most promi¬ 
nently from the electrostatic term 2'nd/Ae in relation 
(25), which is entirely absent in Refs. mm- The more 
complete treatment of our work yields a compressibility 
that is positive for the parameter set of the LAO/STO 
heterostructure and therefore retains the stability of the 
electronic interface-surface system. 


There is also a more technical difference: here, the 
compressibility is found from the second derivative of 
the free energy (see relation (§). In Refs. mm, how¬ 
ever, the compressibility is calculated from the deriva¬ 
tive of the chemical potential, which was identified from 
a model with fixed, density-independent RSOC. Only for 
the derivative of the chemical potential, which yields 
the RSOC was generalized to include the density depen¬ 
dence. Evidently, the two approaches lead to diverse re¬ 
sults even when the analysis is restricted to a single layer 
model without surface term (see Appendix j^. 

In this article we addressed explicitly the dependence 
of the capacitance of the heterostructure on the RSOC. 
As reiterated above, the differential capacitance Cdis is 
not the compressibility in this scheme. The differential 
capacitance is technically found from the second deriva¬ 
tive of the free energy with respect to the charge density 
n at the interface. So Cdiff is an internal capacitance 
of the heterostructure, the capacitance between interface 
and surface. A 3D conducting film may be deposited on 
the surface to have a metallic electrode (gold or YBCO, 
see Appendix]^. The kinetic contribution to the capac¬ 
itance suppresses the capacitance below its geometrical 
value when exchange effects and RSOC can be neglected. 
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FIG. 9. (Color online) Compressibility k (a), interface electronic density n (b), and capacitance Cdiff (c), each in dependence 
on the polar charge density Q. Here me = 0.7me, m^ = 15 me, mi = me, d/e = 10 as, Aso = lOmeV, Ae = 50meV and 
Az.o = 10“^ e^. The parameter Az,i (for the hybridization between d^y- and dxz,yz-orbitals) controls the density dependence 
of the RSOC. The voltage Vo allows to tune the interface charge density. In order to attain a charge density compatible with 
experiments one has to choose a negative voltage; the lowest value of Vo for given Az,i corresponds to an interface electronic 
density of 2 • 10^^/cm^; it is nearly identical for the three sets. The value of Q = le“/uc for the LAO polar film is specially 
marked in the plots. At the point where the system undergoes a transition from n = 0to0<n<Q, the inverse compressibility 
displays a sharp drop as a function of Q. Note that the differential capacitance is only defined for 0 < n < Q, that is, when a 
gate bias can transfer charges from the surface to the interface. 


However, we find that a strong RSOC can enhance the 
capacitance with increasing gate bias to values beyond 
the geometric capacitance. Hence, experiments can ver¬ 
ify a RSOC through a capacitance enhancement. 


A 6-band model is more realistic to investigate the 
capacitance of LAO/STO heterostructures (Sec. IV). 
Most importantly, it allows to derive a Rashba-like SOC 
through the coupling of an atomic SOC (at the Ti sites) 
to the t 2 g bands in the presence of field controlled in¬ 
version symmetry breaking at the interface. The 2-band 
model is interpreted as an effective model for sufficiently 
small filling when only one t 2 g orbital is occupied and the 
splitting of the bands is still linear in momentum k, that 
is, “Rashba-like”. However, already for density indepen¬ 
dent RSOC and low filling, the capacitance as a function 
of bias voltage V behaves differently than in the 2-band 
model. The capacitance increases perceivably (Fig.j^b)), 
dashed brown line) for the 6-band model whereas it is 
constant with respect to V for the 2-band model in the 
absence of a density dependence of the RSOC. In order to 
find capacitances larger than the geometric capacitance 


a 


geom5 


one has to include the density dependence of the 


band structure. 


However, one has to realize that in the regime of 
stronger gate bias V, a variable number of bands are oc¬ 
cupied (see Fig. I^a)) and the splitting of the SOC band 
pairs is not anymore linear in k (see Fig. the middle 


band pair has anyway no linear splitting, even for k ^ 0). 
Then the capacitance as a function of V depends strongly 
on the chosen parameter set (t 2 g band splitting, effective 
masses, SOC, orbital deformation). We could identify 
bias-induced capacitance enhancements of the order of 
5% for parameter sets compatible with the LAO/STO 
heterostructure, but a 30% enhancement as observed in 
experiments [1] appears to be beyond the present mod¬ 
elling. It should be noted that those experiments are in 
a regime of very low charge carrier densities and, there¬ 
fore, it is expected that electronic exchange is much more 
prominent, an effect which would have to be included. 

For a realistic value of the effective distance between 
interface and top gate, d/e = 10 ob, the total compress¬ 
ibility of the system stays positive even for large values of 
the parameter that controls the density dependence 
of the hybridization between the dxy orbital and the dxz.yz 
orbitals. There is no instability which produces a phase 
separation on account of a Rashba-like SOC. 


We portray the basic capacative properties of a polar 
heterostructure in Fig. by relating the compressibil¬ 
ity, the relative value of the interface electronic densitiy 
and the differential capacitance to a variable polar charge 
density. The value oi Q = le“/uc for the LAO/STO het¬ 
erostructure is marked. We assigned Az_i = 33e^a| to 
the largest value of the RSOC that is still compatible with 
experiments (see Sec. IVI. The potential Vb controls the 
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shift of the interface t 2 g-bands with respect to the surface 
states: the lowest value of Vb (for each of the three sets 
of curves for given in Fig. adjusts the value of 
the interface electronic density n to 2 x 10 ^^/cm^, which 
is in the physical range of electronic densities. This cor¬ 
responds to n/Q ~ 0.03 in Fig. [^b). It is obvious from 
Fig.|^a) that for Q = le“/uc and for the LAO/STO val¬ 
ues of Vo I tbe compressibility is well in a positive range. 
However, it is conceivable that for other heterostructures 
with a polarity of Q = 2e“/uc and large RSOC, that is 
sizeable i, a phase separation is attainable with this 
mechanism and a considerable enhancement of the dif¬ 
ferential capacitance over the geometrical value can be 
achieved (cf. Fig. 0 c)). 


ACKNOWLEDGEMENTS 


This work was supported by the Deutsche Forschungs- 
gemeinschaft through TRR 80. We are grateful to J. 
Mannhart for his continuing support of our research. The 
authors acknowledge helpful discussions with R. Ashoori, 
S. Graser, M. Grilli, J. Junquera, A. P. Kampf, Lu Li, 
and V. Tinkl. 


Appendix A: Calculation of the Energies 


1. Single Layers 


For a single layer the strength of the RSOC or depends 
on external fields only, so that or is not a function of 
the density n in that layer. Hence we obtain the relation 
between number of electrons N and chemical potential 
via the energy dispersion, Eq. Q: 


N= 0(M-c^fc)) 

k,A=± 


A 

27r 



dk-kj^e 

A=± 


M - 


f (fc - AttAoPr)^ 
27tNo 



With the abbreviation ( = —2e“;„ = 7 rAoaR /2 and as¬ 
suming fi + C > 0 (otherwise the Heaviside functions and 
hence the integrals are zero), we get 


N = 


A 

27r 



dk-kj^e 

A=± 


+ (7rAoaR)2 


-\k- AttAgOrI 


(A2) 


In our analysis we have to distinguish between fi > 0 
and /X < 0 , where both bands or one band is occupied, 
respectively: 


^ < 0 : 


47rn = 


|•^TNoaR + ^/2^rNofJ,+{^TNoaR)^ 
/ ^rNoaR — ^y2^rNo^t+{^TNoaR)^ 


dk ■ 2k 


^ ^l = 


47rAoQ:R\/27rAo^ -I- 

V? fa 2 

2 ^^^ - 2 = J2N^ 


C 


(A3a) 

(A3b) 


^ > 0 : 


dnn = 


f-s^2TTNofj.+{TrNoaR)^+'wNoaR 


dk ■ 2k 


ff^2ifNotJ.+{TTNoaR,)^-TrNoaR, 


+ 


dk ■ 2k 


= 47rAo/i -I- 4 (ttAgOr)^ 
n 9 n 


(A3c) 

(A3d) 


The same distinction has to be made for the calculation 
of the free energy. 

^ < 0 : 


F<= Y. eHk)e{fi-e\k)) 

k,A=± 

^ nTTNoaR,+f/2TfNofJ.+ iTrNocyRy 


dk 


2 '^ J'KNoa-R — ^2-KNo^+{TTNoaR)‘^ y 27r A( 




— ttRfc' 


2 A _ 

= - TT^Ag Or) V27rAo/7+77rAoQ(R^ 

(A4a) 
(A4b) 


= A 




V6^A3 


O'^R 




Note that t he A = — 1 term is zero, since the 

square root •\/27rAo^ -I- (TrAgaR)^ is always smaller than 
\k + ttAocarI for negative /x. For ^ > 0 this term yields 
the integral with the smaller upper bound: 

^ > 0 : 




A rfy2'7TNotJ,+ i'^NoaRy+TTNoaR. 


dk 


+ 




27r 


dfc 


2-kNo 

k^ 

2nNr 


^ ^ “ xrAoa^n -f y A^al^ 


— ttRfc^ 
+ ClRfc^ 

(A5a) 

(A5b) 


Appendix B: Electrostatic Energy 

We calculate the electrostatic energy of a system sim¬ 
ilar to a LAO/STO multilayer, cf. Fig. before 
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electronic reconstruction positive and negative charged 
planes alternate in z-direction and the lowest layer is 
both electrical neutral and free of electrical fields. If an 
electron charge density —en is transferred from the top 
to the bottom layer (Fig. [^a)), an additional field pro¬ 
portional to n is superimposed. Note that for the RSOC 
only the asymmetric fields are relevant, so that the field 
of the electrons in the bottom layer does not contribute 
to RSOC there. Even if there is only a field in the up¬ 
per half-space, the field important for SOC in the plane 
is the field originating from the surface in the complete 
space. 

The electrostatic energy of a system with I alternating 
(double) layers can be calculated via the energy stored 
in the fields: the electric field if of a planar charge density 
ecr (where e is the electronic charge) is given by 


27r 

E = —ecr 
e 


(Bl) 


and the energy per area between two plates with charge 
density Eea and a distance S between them is 


FesM = 


1 e / drr 

2 An 


-ecr 


= 


2ne'^d} 


a\ 


(B2) 


Hence the energy per area of I double layers with charge 
density ±eQ in an external field, which is generated by 
a charge density ±en, is 


iz /^ _ r/o ^2 , 27re2(l4) 2 

tes/A — - \Q — n) -\ -n 

Si £2 

Fig-ib) shows an equivalent circuit for 1 = 3. 


(B3) 


Appendix C: Three-dimensional layer on Top of 
RSOC Layer 


In some experimental setups, a three-dimensional (3D) 
top gate is deposited on top of the LAO/STO het¬ 
erostructure (see, for example. Ref. my The electrode 
connects with the topmost AIO 2 layer and the electrons 
therein should be treated as a 3D electron system. 

This setup requires to consider an inhomogeneous sys¬ 
tem with field penetration in the top electrode as in 
Ref. [24l We assume a space independent DOS N{Ep) 
within the electrode on the ground that this does not 
modify the qualitative aspects (in fact, modifications are 
tiny if the electronic system is not close to a van Hove 
singularity). Then the kinetic term from the electrode 
free energy Fi translates into a correction of the effective 
distance in the geometric capacitance: the effective dis¬ 
tance d/s is increased by Atf, the screening length (see 
Refs. Ha EH [Ml IMS]). So instead of implementing an 
inhomogeneous electron system in the energy F) of the 
top electrode, one may cast the electrostatic energy into 
the form (cf. Eq. (§): 


Fes/A = 


nde^ 


- + 2Atf ) e" (Q - ny . (Cl) 



FIG. 10. (Color online) Differential capacitance per area 
(solid lines) for a gold electrode on top of the LAO/STO het¬ 
erostructure which is terminated by an AIO 2 layer. Dashed 
lines refer to C = ^ 7 ^. Here msTO = 0.7me, d/e = lOan, 
Q = le“/uc, and oq = 0. Vo ensures that the density on the 
interface is 2- 10^^cm“^ (see Eq. ^ for the significance of Vo). 
The gold electrode has a thickness of ISOan, ruAu ~ 1.1 me, 
Atf = 1.57 an and the capacitance Ca=o/A = 6.67 • 10“® 
for zero RSOC. The system is in phase P2 for all plotted val¬ 
ues of «!, except for the Oi = 10 e^a| curve, where the system 
is in phase Pl'^ to the right of the kink at V = S.Se/as. The 
horizontal dashed brown line indicates the geometric capaci¬ 
tance. 


We assume here that the number of electrons in the 3D 
film is much larger than the number of unoccupied sites 
in the interface layer with RSOC, so that even though 
charge is transferred to the bottom layer, the volume den¬ 
sity nv in the metal remains approximately unchanged. 
This condition can always be fulfilled by a sufficiently 
thick metallic layer. Then the screening length is approx¬ 
imately constant—characteristic of the used electrode 
material. 

In Thomas-Eermi theory Atf depends on the volume 
electron density uy and the dielectric constant Sc of the 
ion cores in the metal: 

^ \lATremiEp) ^ \/ 6ne^nv ’ 

where N{Ep) is the 3D DOS at the Fermi energy. Table|T] 
shows the screening lengths for Au, Ag and YBa 2 Cu 307 . 

Note that Q and ny are fixed, as they are determined 
by the polarity of the LAO layers and the choice of the 
top gate film. In contrast, the electron density n is found 
from the minimization of the free energy. 

In Fig. the differential capacitance Cdia and C = 
eAn/AE are displayed for Au as top gate. We assumed 


e 
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a screening length of I.STob for gold. For the weakest where we suppress the label for the variable Q in Cdiff, 


RSOC, ai = O-le^Og, the capacitance is not altered no- or from the Legendre transform Eq. (D2|: 


tably in comparison to the case without RSOC. For the 


largest value oi = 10e^a| (see estimation in Eq. (16l), 


we find a capacitance increase of approximately 8 % in a 
voltage range 0-100 V. This is the capacitance increase 
with respect to the case of zero RSOC. The increase with 
respect to the geometrical capacitance is 5 % for a gold 
top electrode. 

Different screening lengths cause sizeable differences in 
magnitude of the capacitances: the screening length in a 
YBa 2 Cu 307 film is ^ 10 ob, and the capacitance increase 
is approximately 3 % in the voltage range 0-100 V. The 
increase with respect to the geometrical capacitance is 
only 1 % for the YBCO top electrode. 


Appendix D: Relation between compressibility and 
capacitance and stability criterion 

In this section we assume that the electron density at 
the interface, n, does not take one of the boundary values, 
0 or Q- Moreover, we assume that the appropriate free 
energy is minimal with respect to n(Q, V) (see below). 

The Helmholtz free energy Etot of the system is a func¬ 
tion of the total electron density Q and the electron den¬ 
sity on the interface n (cf. Eq. (§): 

FtotiQ, n) = Eo(n) + Fi{Q - n) + F^siQ, n) - eVonA 

The partial derivative of FtotiQ,^) with respect to n is 
the voltage between the plates: 


dFt, 


dn 


= eAV {Q, n) 


(Dl) 


This thermodynamic relation may be used to identify 
^(Qi F). We introduce a Legendre transformation of F^ot 
with respect to the variable n: 


Cdiff(E)M = 


e dn{Q, V) 
dV 


1 d'^F 
'~AdV^ 


(D4) 


Of course, the two capacitances from Eqs. (D3) and (D4) 
are identical in the sense that Cdiff(n(E)) = Cdiff(E) 
holds. 

On the other hand, the inverse compressibility 
V) in dependence on an external voltage V is the 
second total derivative of F with respect to the total 
electron density: 


-1 /i/zr2 d2F(Q,n(Q,E)) 


■ w = 


dQ 2 


(D5a) 


= dlF + 2dr,dQF {dgn) + dlF {dgnf 

+ dnF {din) (D5b) 

Note that dnF = dF/dn is not a derivative with respect 
to a thermodynamic variable. In fact, the thermody¬ 
namic variables are Q and V for F. The partial deriva¬ 
tive of F with respect to n is identified from the rhs of 
Eq. (D2). The thermodynamic potential F is minimal 
with respect to all values of the internal variable n, a 
necessary condition for thermodynamic stability. 

According to Eqs. (D2) and (Dl), the last term on the 


right hand side of Eq. (D5a) is zero. Under a further as¬ 


sumption we can derive more elaborate relations between 
differential capacitance and compressibility: the depen¬ 
dence of the total free energy Ftot on the total electron 
density Q shall be exclusively of the form Q — n, so that 
also 


d^F 

W 


d'^F 

dQdn 


(D6) 


F{Q, V) = Ftot(0, n(Q, V)) - eA U • n{Q, U) (D2) 

where F{Q,V) depends on the total electron density Q 
and the external potential V between interface and sur¬ 
face electrode. 

The differential capacitance can be derived either from 
Eq. (pi) 


A/C'diff(n) = 


dV{Q,n) 
e dn 


1 d^Fto, 

A dn^ 


(D3) 


TABLE I. Material properties and resulting screening lengths 
Atf for gold, silver and YBCO. Here, nv is the volume elec¬ 
tron density and £c the dielectric constant of the ion cores in 
the metal. 



Ef [eV] ec 

nv [ 10 ^^ cm 

Atf 

An 

5.53 [351136] 6.9 [3l| 

59.0 [351136] 

1.57 

Ag 

5.49 [SnilSS] 2.5 [31] 

58.6 |35l I36| 

2.99 

YBa2Cu307 

1.0 [38] 25 HHj 

2. ..6 [38] 

7.4. ..12.8 


is valid. This assumption is reasonable, since Q — n is 
the electron density in the surface layer, which should be 
the only density dependent quantity that determines the 
energy of the surface system. Hence we only demand the 
electrostatic energy not to have some unusual form. Note 
that, while Etot(Q, n) = Fto tjQ — n, n) is sufficient for the 
derivation below, only Eq. ( |D6[ ) is necessary. Under the 
validness of this relation the compressibility becomes: 




d^F 

dQdn \ dQ ) 


d^F / dn\ 

dn^ \ dQ ) 

(D7) 


We use the total derivative of the first condition Eq. (|D1|) , 


d dF 
dQ dn 

d^F d^F / \ 

dQdn dn^ \ dQ J ’ 


(D8) 
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to replace dqdnF in the compressibility above: 


-1 


= S (^ 


dn 
drF \dQ 

dn 

C'diff dQ 


dn 

~ ^ 


dn 

~ ^ 


(D9) 


This relation between the compressibility and differential 
capacitance can be further analyzed to deduce a stabil¬ 
ity criterion, i.e., the condition for positive compressibil¬ 
ity. Sinc e we assumed that the density n obtained from 
Eq. (Dl) yields a minimum of the free energy, 


d'^F 

dn^ 


> 0 , 


(DIO) 


the differential capacitance C'diff is positive, which is a 
necessary criterion for the stability. Note that the second 
partial derivatives of F and Clot (with respect to n and 
Q) are equal. Hence there is an addtional condition to 
keep the compressibility positive: 


At > 0 




(Dll) 


The system is unstable if adding of electrons causes ei¬ 
ther electrons to flow from the interface to the surface 
{dqn < 0) or more than the added electrons to flow to 
the interface {dqn >1). As we will show below, the 
first alternative does not occur in our model, dqn < 0 is 
in principle the reverse case of dqn > 1, i.e., when the 
isolated surface system would display a negative com¬ 
pressibility. 


Furthermore, Eq. (D8| can be solved for dqn, 


dn _ d^F /d'^F 
dQ dQdn / dn'^ ' 


(D12) 


This result inserted into the compressibility Eq. (D9) 
yields: 


A/Q^ = - 


d^F 
dQdn 
d^F 

W 

d^F 

W 


d^F 


1 - 


dQdn 
d^F /d^F\ 
) 


d‘^F \ 
dn^ J 


( C'diff 

e2A2 dQ^ J 


(D13a) 

(D13b) 

(D13c) 


We consider a system where the electrons in the sur¬ 
face electrode are not subject to a RSOC or other effects 
which may produce a negative compressibility. Corre¬ 
spondingly, we may assume that d'^FijdQ'^ > 0. As the 
electrostatic contribution to the free energy also gener¬ 
ates a positive second derivative with respect to Q, the 
free energy F has to obey 


d^F 


> 0 . 


(DM) 


Since dndqF is negative (Eq. (D6) with relation (D14)) 
and d‘^F /dn^ is positive in the considered regime, 
Eq. ( |D8[ ) shows that dqn has to be positive in our model. 

We consider the regime of positive capacitance—is it 
then possible to identify a negative compressibility of 
the comple te electronic system? In fact, according to 
Eq. (D13b) and the positivity of d’^FjdQ’^, the compress¬ 


ibility is negative if 


K < 0 


d'^F ^ d'^F 


dn^ 

d\Fo 


< 


d^Ftoi 

dQ^ 

Fi)/A 


. irde^ 


dn^ 

d^Fo/A 

dn^ 


■nde^ 

2- < 0 


d'^FxjA „7rc?e^ 

dn^ e 

(D15) 


This inequality does not involve properties of the surface 
system explicitly. With the assumption that C^is is posi¬ 
tive, one finds that the compressibility of the heterostruc¬ 
ture is negative, if the sum of the inverse compressibility 
of the interface system (times n^) and half the inverse 
geometric capacitance 2'!Tde^/e is smaller than zero. 

We emphasize that the inverse capacitance is a sum, 
the terms of which are generated by the second deriva¬ 
tives of different contributions to the Helmholtz free en¬ 


ergy Etot in Eq. (D3), consistent with previous analy¬ 
ses (e.g.. Ref. mr The inverse compressibility is not 
the sum of these term: it is the asymmetric setup that 
cause s (Cdiff/A)”^ ^ K“^/(e^Q^). However, the relation 
(D15) contradicts the identification of a negative com¬ 


pressibility in Refs. mm- there the authors missed the 
“electrostatic compressibility” 27rc?e^/e entirely and the 
compressibilty of the interface electronic system was eval¬ 
uated differently (see Appendix]^. 

In our model the total free energy Ftot is of the form 
presented in Eq. (|5l), which reads with the abbreviation 


D = nde'^/e = (eQA) ■ A/C^eom- 

Ftot{Q -n,n)= Fo(n) -H {Q - nf A 

+ Dn^A - eVonA. (D16) 

Due to the analytical form of dqFtot = (l/A^i + ‘^F>)A, 


differentia 
1 


Eq. (DI3b) yields a direct relation between the measured 
capacitance and the total compressibility: 


K ^ = 


1 


2D 


1 - 


C'diff f 1 


e’^A \Ni 


2D 


and the transfer of charge to the interface is: 
dn _ / 1 1 A \ Cdiff 

dQ~ WNx 2Cg,,^) A • 


(D17) 


(D18) 


Obviously, for the symmetric case without RSOC and 
with equal DOS at the interface and surface electrode. 
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FIG. 11. (Color online) Differential charge transfer dn/dQ 
(a), and the corresponding interface electronic density n (b), 
both with respect to Az,i, which tunes the density dependence 
of the RSOC. The density n is controllel by Vq. Here, me,= 
0.7nie, = 15me, mi = me, d/e = 10ob, Q = le“/uc = 
6.6 • 10^'‘/cm^, Aso = lOmeV, Ae = 50meV, and Az,o = 
10“^ e^. We find Cdiff > Cgeom for dQti > 41/80. The jump 
to the horizontal line with dn/dQ = 1 in (a) takes place when 
all electronic charge is accumulated at the interface, see (b). 


No = Ni, one finds dn/dQ = | and, consequently, n^K = 

Cdiff/ 

For effective interface-surface distances d/e > 10 ob, 
the electrostatic term 2D = ^ • A/Cgeom dominates over 
the inverse compressiblities l/7Vo,i, and the charge trans¬ 
fer dn/dQ is close to ^ for vanishing RSOC. Therefore, 
sizeable deviations of dn/dQ from ^ can be indicative 
of a substantial RSOC in these systems. With the mea¬ 
surement of Cdiff and knowledge of Cgeo m and Ni o ne can 
determine dn/dQ and k through Eqs. (D18) and (D17l, 
respectively. 


The electronic transfer dn/dQ is displayed in Fig. 11 


as 


a function of the density dependent hybrization, which ef¬ 
fectively controls the strength of the RSOC (see Sec. IV). 
The limiting value of dn/dQ = 1 is reached, when for 
large RSOC all additional charge carriers are accumu¬ 
lated at the interface (horizontal lines in Fig. 11). 


Appendix E: Compressibility for the phase with 
charge-depleted surface (Q = n) 


ity We take the strength of the RSOC or, density 
dependent, and the density dependence is implemented 
already in the free energy. However, Caprara et al. [5] de¬ 
termine the chemical potential from a model with given 
Or, then introduce the density dependence of or at that 
level and calculate the inverse compressibility from 
the first derivate of the chemical potential with respect 
to n. 

Here we compare the resulting compressibilities n and 
Kc of these two approaches. Therefor we use the abbre¬ 
viation C = 7rVoQ;^/2, which equals the lowest energy in 
the dispersion Eq. ([^, and rewrite the energy Eq. ^ of 
a layer with RSOC: 

f n^ 

(El) 

Eor density independent or and hence density indepen¬ 
dent Q the chemical potential /i is given by 




= = ) 4:N^C 


- C < 0 


Adn 


Nr 


-2C ,M>0 


(E2) 


Caprara et al. [5] then set ckr = Q!R(n), and their inverse 
compressibility is 


( 4 ) 


-1 


= ri^^=n4 i^Noc) 


= n = n 
dn 


- 2 C' 

No ^ 


(E3) 


Here Q = dC,/dn and C,” = d’^(//dn^. In contrast, if we 


set Or = aR(n) already in the energy in Eq. (El): 


= 


dFjf 

Adn 



,fi<0 

,H>0 


. (E4) 


(«:<)-^ = («:<)-! + 




n n 


C' <' 


\2N^C,^ V 3C 2 6 

- iC - C'n) (E5a) 


(k^) ^ — (kq) ^ -I- (—^C + 7;No ((CO^ + CC") ~ 

(E5b) 


Eor aR = ain this simplifies to 


The first derivative of the energy with respect to the 
electron density n yields the chemical potential while 
the second derivative generates the inverse compressibil- 


1 ( -TrNoaln ,/r < 0 

(k>) (E6) 
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and 



1 

iVo 


—SttA^oQ^i’^ I M < 0 

(!>TTNi)a\n + 271^ ,/i>0 

(E7) 


As one can see from Eqs. (E5b) and (E7|, (kq)“^ and 
differ by a factor of 3. The deviation for 
is less since I/Nq dominates (E5b) and (E7) for realistic 
values of n, ai, and Nq. 
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